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We investigate the unpolarized pion and kaon fragmentation functions using the nonlocal chiral- 
quark model. In this model the interactions between the quarks and pseudoscalar mesons is man- 
ifested nonlocally. In addition, the explicit flavor SU(3) symmetry breaking effect is taken into 
account in terms of the current quark masses. The results of our model are evaluated to higher 
Q^ value Q^ — 4GeV'^ by the DGLAP evolution. Then we compare them with the empirical 
parametrizations. We find that our results are in relatively good agreement with the empirical 
^T^ parametrizations and the other theoretical estimations. 
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I. INTRODUCTION 



Fragmentation functions are important ingredient for understanding the structure of the hadrons, because they play 
T^ a crucial role in analyzing the processes involving hadrons. For example, one needs the unpolarized fragmentation 
Mh functions to analyze the semi-inclusive processes in the electron-positron scattering, deep-inelastic proton-proton 
~^ scattering, and so on [TMllj. Furthermore, to extract the chiral-odd transversity parton distribution of the nucleon 
one needs more complicated fragmentation functions, such as the polarized dihadron fragmentation functions and 
the Collins fragmentation functions. Because of their fundamental importance, those functions have been studied 
intensively for decades but still not fully understood yet. 

The unpolarized fragmentation function D^{z) represents the probability for a quark q to emit a hadron H with 
Xj the light-cone momentum fraction z. It can be written with the light-cone coordinate as follows [l2l[T3] : 

in 1 I- 

p Dfiz,k^^,^i) = -J dk+Tr[Aik,p,tih-]U-^p-. (1) 

^rr Here, k± — {ko ± ks)/^^ and the correlation A(fc,p, fi) is defined as 

^ Aik,p,f,)=J2 f ^^e+^''-HomO\H,X){H,XmO)\0), (2) 

^ where k, p indicate the four- momenta for the initial quark and fragmented hadron, respectively. Furthermore, z is the 
^ longitudinal light-cone momentum fraction possessed by the hadron and ^ denotes a renormalization scale at which 
the fragmentation process is computed. All the calculations are carried out in the frame where the z-axis is chosen 
to be the direction of fe. Consequently the transverse momentum of the initial quark fej^ is zero in this frame. On 
the other hand, kx = k — [{k ■ p)/|pp] p, defined as the transverse momentum of the initial quark with respect to the 
direction of the momentum of the produced hadron, is nonzero. The integrated fragmentation function satisfies the 
momentum sum rule: 



/ y2zD^{z,^i)dz = l, (3) 
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FIG. 1: Quark fragmentation cascade process. 

where H represents all the fragmented hadrons. Eq. (|3| means that all of the momentum of the initial quark q is 
transferred into the momenta of the fragmented hadrons. Empirically, information of D^{z) has to be extracted 
from the available high-energy lepton-scattering data by global analysis with appropriate parameterizations satisfying 
certain constraints [TH - fT^ . 

From theoretical points of view, it is impossible to study fragmentation functions directly by lattice QCD because 
they are defined in Minkowski space. Nevertheless, there have been numerous works for the fragmentation functions 
based on the effective QCD models so far. In Ref. pO], the Nambu-Jona-Lasinio (NJL) model has been used to 
calculate the fragmentation functions. Monte-Carlo (MC) simulations with supersymmetric QCD were also carried 
out to obtain the fragmentation function up to a very high energy in the center-of-mass frame y/s [3T]. The Collins 
fragmentation functions, which play an essential role in the transverse-spin physics, have been studied as well in the 
quark-pseudoscalar (PS) meson coupling model [121 HSl [22] . Note that dihadron fragmentation functions have been 
investigated in the same theoretical formalism |23j . For brevity, we will simply call the unpolarizcd fragmentation 
functions as the fragmentation functions from now on. 

In Rcfs. |24l[25], we have already employed the nonlocal chiral quark model (NLChQM) with the explicit flavor 
SU(3) symmetry breaking to calculate the elementary fragmentation functions. This instanton-motivated approaches 
were used for computing the quark distribution amplitudes, manifesting the nonlocal quark-pseudoscalar (PS) meson 
interactions P5ti29j . NLChQM have been applied to determine various nonperturbative quantities and obtained 
the results which are in good agreement with experiments as well as lattice QCD (LQCD) simulations [30ti35). 
The elementary fragmentation functions calculated in Refs. [23J [5S] are the functions in Eq.S with the following 
approximation: 

J2 \KX){h,X\ « ^ |/i = gQ, X = Q){h = qQ,X = Q\. (4) 

X h 

Here h denotes the PS meson. In other words, we just calculate the one-step fragmentation process: q ^ h = {qQ) + Q- 
In Refs. j24l I24j . the renormalizcd fragmentation functions are obtained by re-scaling the elementary fragmentation 
functions. After DGLAP evolution the results of the renormalizcd fragmentation functions agree with the empirical 
ones reasonably well except at small z. However such renormalizcd fragmentation functions do not satisfy the sum 
rules of Eq. ([3|). Furthermore some of the fragmentation functions , such as u — > tt~ and s —>■ K~^, are zero because 
their elementary fragmentation functions are identical zero. Hence it is necessary to include the quark-jet contribution 
from many-step fragmentation processes. 

In Refs. pm [3BHSH]j the NJL model has been applied for the fragmentation functions including the quark-jets 
and resonances. The momentum sum rules is satisfied automatically according to their approach. It turned out 
that the quark-jet contributions provide considerable contributions to the various fragmentation functions at the 
small z region. The approach in [501 [55H55] is actually applicable for any effective model. With this method, one can 



generate the fragmentation functions D^[z) by solving the coupled differential- integral equations where the elementary 
fragmentation functions dh{z) appear in the kernel. In this article, we will extend our previous works in |241 [25j by 
including the quark-jet contribution of the fragmentation functions using the method developed in pUl [5S1 - I55] . We 
will evolved our results to higher Q^ values. Then we will compare our results with the empirical parameterizations 
and the results from the NJL-jet model. 

The present work is organized as follows: In Section II, we briefly sketch our results of the elementary fragmentation 
functions based on NLChQM. In Section HI, we explain the method of calculating the quark-jet contribution to the 
fragmentation functions and apply this method to our model. In Section IV, we present and discuss our numerical 
results of fragmentation functions at Q^ — 4GeV^. The final Section is devoted for the conclusions and future 
perspectives. 

II. ELEMENTARY FRAGMENTATION FUNCTIONS IN NLCHQM 

In this section, we briefly explain how to derive the elementary fragmentation functions in NLChQM. This model is 
motivated from the dilute instanton liquid model(DILM) PUH15] in which the nonperturbative QCD effects are from 



the nontrivial quark-instanton interactions in the dilute instanton ensemble. However DILM is defined in Euclidean 
space since the (anti)instantons are well defined there as the tunneling between the infinitely degenerate QCD vacua. 
However, there are still some works which generalize DILM to calculate some physical quantities such as the light-cone 
wave function which are properly defined only in Minkowski space [26 29 . Following those studies, we adopt the 
efi'ective chiral action (EChA) from NLChQM in Minkowski space as follows: 



Scs[mq,h] = -iSpln 



where Sp and rhg denote the functional trace Tr J d* 



m/ - \ImX^)W'' \Im{^) 



(5) 



over all the relevant spin spaces and SU(3) current- 



quark mass matrix, diag(m„,r?T,d,''7is), respectively. Here we choose the following values: (m„,md,TOs) = (5,5,150) 
MeV. As mentioned in Refs. PH |2S], to derive EChA Eq. ^ we simply replace the Euclidean metric for the 
(anti) instanton effective chiral action by the one of Minkowski space. The interactions between the quarks and 
the nonperturbative QCD vacuum generate the momentum-dependent effective quark mass which can be written in 
a simple n-pole type form factor as follows 



M{d'^) = Mo 



iA2 - 92 + ie 



(6) 



where n indicates a positive integer number. We will choose n = 2 as in the instanton model [28l|29]. In Eq. ^ 
A « /i stands for the model renormalization scale. It is related to the average (anti)instanton size p in principle and 
takes the value A « 600 MeV |13]. The nonlinear PS-meson field, i.e. U^'^ takes a simple form (i.e. Ref. [33]) with the 
normalization chosen to be consistent with the definition of the fragmentation function in Eq. (IT]) [T^] : 



W' [h) = exp 



»75(A-fc) 
2Fh 



= 1 



i75(A./i) {\-hf 



2Fh 



8Fi 



(7) 



where Fh and A" stand for the weak-decay constant for the PS meson h and the Gell-Mann matrix. The flavor SU(3) 
octet PS-meson fields are given as: 



Xh = V2 



7I-° 



■/6 



r] 






RO 



_V6' 




(8) 



By expanding the nonlinear PS-meson field from EChA in Eq. ([5]), one derives the following effective interaction 
Lagrangian density in the coordinate space for the nonlocal quark-quark-PS meson vertex: 



-qQh 



2Fh 



Q\^M{d^)j,{X-h)^jM{d^)q 



(9) 



As a result, we reach a concise expression for the elementary fragmentation function q{k) — > h{p) + Q{r) from 
NLChQM: 



d'Hz,kl,^i) = 



c: 



MkMr 



87r3;2(l-z) 2Fl 



z{P - Ml) + (fc2 + Ml - 2MqMQ -2k-p) 

(P - M2)2 



(10) 



where C^ indicates the flavor factor for the corresponding fragmentation processes and is given in Table u\ in the 
Appendix. Notice the flavor for the initial quark q in Eq. (10) is written explicitly. Furthermore, the momentum 
dependent effective quark mass reads: 



Mf = M 



2A2 



n 2 



2A2 - £2 



(11) 



Here we have used the notation: Mq = mq + Mo . The value of Mq can be fixed self-consistently within the instanton 
model [5D1 - I551 [551 - H5] with the phenomenological (anti)instanton parameters p ~ 1/3 fm and _R w 1 fm. This will lead 
to Mq « 350 MeV. The masses for the pion and kaon are chosen to be TOtt.x = (140, 495) MeV throughout the present 
work. Taking all the considerations into account, one arrives at a concise expression for the elementary fragmentation 
functions: 
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, krp, P) 



Q MkMr 

87r3 2Fi 



[z^kl + [{z - 



l)Mq + MqX' 



[z^kl + z{z - 1)M? + zMl + [l- z)mlf '' 



(12) 
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FIG. 2: The fragmentation functions of zDZ, {z) (upper panel, left), zDZ, {z) (upper panel, right), zD^ {z) (bottom panel, 
left) and zD^ (z) (bottom panel, right). The dashed lines denote the results of NJL model. The solid lines stand for the 
results of nonlocal chiral quark model employed in our calculations. 



where Mk and Mr are the momentum-dependent quark masses manifesting the nonlocal quark-PS meson interactions: 

AU2K^z{l-z)f ,., Afo(2A2)2 



Mk 



[z^kl + z{z - 1)(2A2 - ,52) + zM^ + (1 - z)mlf ' 



Mr 



(2A2-A?2)2- 



(13) 



As in Ref. |24j . a free and finite- valued parameter b has been introduced in the denominator to avoid the unphysical 
singularities. Notice the singularities arise in the vicinity of (z, A;t) = due the present parametrization of the effective 
quark mass as in Eq. (ro|. At the renormalization scale in our model, the elementary fragmentation function can be 
evaluated further by integrating Eq. (12 1 over /ct: 



<(z,/i) 



/•oo 

27rz^ / d'^{z,k'^,fj,)kTdkT. 
Jo 



(14) 



III. FRAGMENTATION FUNCTIONS WITH THE QUARK-JET CONTRIBUTION 

To calculate the quark-jet contribution to the fragmentation functions within our model, we follow the approach in 
Refs. [20 1I36H38] . The elementary fragmentation functions dn{z) are re-defined as follows, 



E/4'w = E/^'?w^^ = 1' 



(15) 



where the complementary fragmentation functions d^{z) are given by 

d^(z) = d^(l-z). h = qi 



(16) 
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FIG. 3: The fragmentation functions of zDZ, (2) (upper panel), zD^ (z)(down panel, left) and zD^ (2) (down panel, right). 
The dashed lines denote the result of NJL model. The solid lines stand for the results of nonlocal chiral quark model employed 
in our calculations. 



The fragmentation functions D^{z) should satisfy the following integral equation: 



Dl{z)dz = dl{z)dz + Y,f dyd'^iy)D'^ (- 



z\ dz 
V ' 



(17) 



Note that D^{z)dz in Eq. ( 17) has an interpretation: D^{z)dz is the probability for a quark q to emit a hadron which 
carries the light-cone momentum fraction from z to z + dz. Since d9{y)dy is the probability for a quark q to emit a 
hadron with flavor composition qQ at one step and the final quark becomes Q with the light-cone momentum fraction 
from y to y + dy. Eq. (17) actually describes a fragmentation cascade process of hadron emissions of a single quark 
depicted in Fig.(fl]). 

To solve the coupled integral equations Eq. (171, we apply the MC method developed in Ref. [37]. Namely, we 
simulate the fragmentation cascade of a quark iVtot times and each time the fragmentation cascade stops after the 
quark emits Minks hadrons. The fragmentation function D^{z) is then extracted through the average number of type 
h hadron with light-cone momentum fraction z to z + Az, N^{z, z + Az), by 



Z?,^z)Az=-^^7V,''(z,z + Az). 

JVtot 



(18) 



When iViinks increases, one finds that £'^(z)Az increases in the low z regime. This is due to the fact that when 
more steps of a fragmentation cascade are considered, more hadrons with low z are emitted. One interesting feature is 
that when z — ^ 0, zd^{z) — >• but zD^{z) — >■ constant. This is true for both of the NJL-jet model and our model. The 
value of £'^(z)Az becomes insensitive to the value of iVtot and Minks when Mot and Minks are large enough, implying 



that the result of the MC simulation converges to the solution of Eq. (17). When Minks reaches 8 the result is already 



convergent in the case of the NJL model [37]. However for our model the results start to converge as Minks > 15. 
This can be explained by the fact that in the NJL model the peaks of the elementary fragmentation functions occur 
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FIG. 4: The fragmentation functions of zD^ (z) (upper panel, left), zD^ {z) (upper panel, right), zD„ [z) (bottom panel, 

left), and zDf {z) (bottom panel, right). The dashed lines denote the results of NJL model. The solid lines stand for the 
results of nonlocal chiral quark model employed in our calculations. 



at higher z value than in NLChQM. It indicates that the probability of a quark carrying mediuni momentum fraction 
(0.4< z <0.8) emitting a hadron is larger in our model than in the NJL model. Hence the MC simulation of NLChQM 
needs the larger value of A^unks to reach convergence. 

In this article the initial quark is only limited to be the light quark, namely, q = u,d, s. In addition, the fragmented 
hadrons only contain pions and kaons. Naively one has 42 fragmentation functions. According to charge conjugation 
and isospin symmetry, there are only 11 independent ones. Notice among these 11 elementary fragmentation functions, 
only four of them are not zero. We call these direct fragmentation functions: 

Dl\z) = Driz) = Df{z) = Dfiz), Dfiz)^Df{z) = Dfiz) = Df{z), 



-iXO 



^RO, 



D^ (z) = D^ iz) = D^ iz) = Df{z), Df {z) = Df" {z) = Df {z) = Df [z). 

In Fig. 2^ we present zD'^ (z), zZ?J (z), zD^ (z) and zDf (z) at /i^ = 0.36 GeV^. The solid lines represent 
the nonlocal chiral quark model (NLChQM) results and the dotted lines stand for the corresponding calculations 
determined from the NJL-jet model. 

Here we briefly discuss the main features of the results of our model and the NJL-jet model. The direct fragmentation 
functions ones are presented in Fig. [2J In the case of u — >■ tt"*" , we find that the shapes of the curves of these two models 
are very different. The peak of our curve occurs at z=0.5, but the peak of the NJL-jet curve takes place at z = 0.8. 
When z > 0.6 our result is smaller than the NJL-jet one. Between z = 0.4 and z = 0.6, NLChQM result increases as 
z decreases but the NJL-jet curves decreases. Below z=0.4 the two curves behave similarly but the magnitude of the 
NLChQM result is about twice of the NJL-jet result. Between z = and z = 0.2 there are plateaus for the both curves. 
In our model, the possibility for the fragmented pion carrying a small momentum fraction (z < 0.4) is much larger 
than in the NJL-jet model. Another direct fragmentation function of pions is Dl^ (z). The elementary fragmentation 
function d^ (z) is exactly one half of d^ (z), dj {z)—^d!^ (z). After including the quark-jet contribution the shapes 
of zD'!^ and zDJ become quite different. Our zDl^ is quite flat between < z < 0.4, then it decreases as z increase 



as z > 0.4. On the contrary, our zD^ is flat between < z < 0.2, then it increases as z increases till z ^ 0.4. 
For z > 0.4 the curve of zDJ^ decreases as z increases. It shows that the eiTect of including the quark-jet is very 
pronounced in our model compared to the NJL-jet model. It is worth mentioning that the relation zDJ « |^« is 
held for z > 0.6 for both of the models. This can be explained as follows: since it is difficult for a quark which has 
already emitted hadrons in the cascade process to be fragmented into a hadron with high momentum fraction, the 
quark-jet contribution is less important at high z regime. For the fragmentation process of s — >■ K~ , we find that 
the quark-jet contribution is small for both of the models. The value of zD^ in the regime z < 0.2 is negligible. 
However, for the fragmentation process of u — >■ K'^, it is clear that our result is much smaller than the NJL-jet result. 
It is of no surprise because our elementary fragmentation function d^ (z) is very small. 

From Table [ll there are several elementary fragmentation functions, such as d^ and dj , are identically zero. After 
including the quark-jet contributions, those fragmentation functions are no longer zero. We call these fragmentation 
functions indirect fragmentation functions and they are generated from the process of fragmentation cascade (i.e. 
Fig. [ij. Those indirect fragmentation functions are listed as follows: 

DZ'iz) = Dfiz) = Df{z)^Dfiz),Driz)^Dfiz)^Dfiz)^Dfiz), 

Df{z) ^ Dr{z)^Dr{z)-Driz),Dr{z)=Dr{z)-Dr(z)^Df{z), 

Df{z) = Df{z) = Df-iz) = Df{z),Df{z) ^ Df{z) = Df{z) = Df{z),Df{z) = of {z). 
For the indirect fragmentation functions of pions depicted in Fig. [31 we first observe that in NLChQM the shape of 

— — — 

zD'Z (z) is somehow similar to zD^ (z). However, the plateau of the zD^ (z) (0 < z < 0.2) is only half of one 
for zD^ (z) (0 < z < 0.4). The magnitudes of zDJ (z) and zDJ (z) are roughly the same at < z < 0.2. As z 

increases zDJ (z) decreases much faster than zD^ does. To compare to the NJL-jet model, we find that in the low 
z regime the NLChQM curve is almost twice larger than the NJL-jet result, but the two curves are very close as 
z > 0.5. The results of zD"^ (z) and zDJ (z) are almost identical. Unlike zD^ , these two fragmentation functions 
monotonically decrease even from very small z. Again the NLChQM results are larger than the NJL-jet results in 
the regime of < z < 0.2. It is because the fragmented 7r+ meson here is enfitted by the multi-step processes, such 
as s — > K^u, u — )• TT+d. Since the peak of df (z) is around z=0.8, the secondary u quark most likely carries small 
momentum fraction z < 0.2. The value of DJ (z = 0.2) of NLChQM is larger than the corresponding one in NJL-jet 
model. As a result the chance of a s quark to be fragmented into a pion is larger in NLChQM. 

In Fig. kjwe present four indirect fragmentation functions of the kaons, zD^ (z), zD^'^ (z), zD^ (z), and zD^ (z). 
The common feature of NLChQM results is that they are all almost one order magnitude smaller than the pion ones. 
In contrast to the pion ones, the NJL-jet results are much larger than ours. It is because these indirect fragmentation 
functions are related to the kaon elementary fragmentation functions which are tiny. Consequently those associated 
indirect kaon fragmentation functions are also suppressed. For example, the process s — >■ K^ is the combination of 
the two processes such as m — > vr+d, d — >■ K^s. Since d^ {z)=d^ (z) is very small, naturally D^ is also very small. 
In summary, we find that our results are substantially different from the NJL-jet model results after including the 
quark-jet contribution. For pion ones the NLChQM results are higher than the NJL-jet ones in the medium and low 
z regime. For the kaon ones the situation is rather different. For s -^ K~ we arrive at a result similar to that from 
the NJL-jet model, but our investigation implies that the process of u -^ K^ is highly suppressed. For the other 
channels NLChQM results are always smaller than the NJL-jet model. 

IV. NUMERICAL RESULTS OF THE FRAGMENTATION FUNCTIONS AT Q^ = 4 GeV^ 

In this section we will present our results at Q^ = 4 GeV^ and compare them with the empirical parametrizations 
and the NJL-jet model resuffs. We employ QCDNUM17 [T71 gSl IMj to evolute our resuhs from Q^ = 0.36 GeV^ 
to Q^ = 4GeV^. Since DJ (z) is the most pronounced process, therefore, the initial momentum for evolution is 
determined by a reasonable agreement between our evolution result of I?J[ (z) with two empirical parameterizations, 
namely the HKNS parametrization [H] and the DSS parametrization [18]. These two empirical parameterizations are 
used for comparison of other fragmentation functions as well. The results of direct fragmentation functions evolved to 
Q^ = A GeV^ are given in Fig. [5J The dashed lines denote the results of NJL model. In addition, the solid lines stand 
for the results of nonlocal chiral quark model employed in our calculations. Finally, the dotted and dot-dashed lines 
are the HKNS curve and DSS curve, respectively. The uncertainty bands are provided by the HKNS parameterization. 
The NLChQM result of zD'^ (z) is within the uncertainty band of HKNS result. In the high z regime our result 
is consistent with the two parameterizations and the NJL-jet model. In the medium z region (0.3 < z < 0.7) the 
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FIG. 5: The fragmentation functions of zDJ (z) (upper panel, right) and zDZ, {z) (upper panel, left), zD^ {z) (bottom panel, 
right) and zDf (z) (bottom panel, left). While the dashed lines denote the results of NJL model, the solid lines stand for the 
results of nonlocal chiral quark model. The dotted lines are the HKNS curves and the dot-dashed lines are DSS curves. HKNS 
and DSS are two empirical parameterizations of the fragmentation functions. The uncertainty bands are according to HKNS 
parameterizations. 



NLChQM result is slight higher than HKNS and DSS. Between z ~ 0.1 and z — 0.4, it appears to be a plateau then 
turns up at z = 0.1. On the contrary, the NL-jet result is slight below the HKNS and DSS parameterizations between 
z = 0.2 and z = 0.4 and turns up at z = 0.2. In the case of zD'^ , NLChQM result is slightly higher than the two 
parameterizations between z=0.3 and z=0.7. On the contrary, the NJL-jet result is clearly below the uncertainty 
between z=0.2 to z=0.4. For the both cases of zD'^ and zD^ , the results of NLChQM and NJL-jet agree with 
HNKS and DSS quite well in the high z regime z > 0.7. We now turn our attention to the case of zD^ . We find 
that both of NLChQM and NJL-jet curves are within the uncertainty band. However, the shapes of the two curves 
are completely different. The NLChQM curve increases from z=0 to z=0.05 then decreases rapidly from 0.05 to 0.2, 
after which it decreases slowly. The NJL-jet curve simply decreases as z increases. Unfortunately, neither of them 
catches the feature of the empirical curves. Model results are too small compared with the empirical ones in the 
region of z > 0.7. Between z=0.4 and 0.7, the NJL-jet model result is comparable with two parameterizations. On 
the contrary, the NLChQM result is always too small. In the low z region, NLChQM curve becomes comparable with 
the empirical ones but the NJL-jet model becomes far too large. The fragmentation function zD^ (z) is probably 
the most problematic one for the model calculations. On the one hand, the NLChQM curve is too large compared 
with the empirical curves in the high z region, but it turns out to lie between the two empirical curves in the region 
of z < 0.4. On the other hand, the NJL-jet curve agrees with DSS curve in the high z region quite well but becomes 
far too large as z < 0.1. Both the model curves fall outside the uncertainty band of HNKS parameterization from 
medium to high z regime. 

Now, let us discuss the indirect fragmentation functions. The ones of pions are depicted in Fig. [6] The most 
successful channel for NLChQM model is zD^ . Between z = 1 and z = 0.2, the NLChQM curve agrees with both 
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FIG. 6: The fragmentation functions of zDZ (z) (upper panel) and zD^ (2:)(down panel, left) and zD^ (z)(down panel right). 
The dashed lines denote the result of NJL model. The solid lines represent the results of nonlocal chiral quark model. The dotted 
lines are the HKNS curves and the dot-dashed lines are the DSS curves. HKNS and DSS are two empirical parameterizations 
of the fragmentation functions. The uncertainty bands are according to the HKNS parameterizations. 

of HKNs and DSS very well. Between z = 0.1 and z = 0.2 the NLChQM curve still coincides with DSS and is a little 
bit higher than HKNS curve. This excellent agreement disppears when z < 0.1. Compared with the NLChQM curve, 
the NJL-jet result locates at the lower edge of the uncertainty band in the high z region and turns upward at z = 0.3. 
For the fragmentation process of s — > tt, both of NLChQM and NJL-jet models give a very similar result. However 
they both underestimate the fragmentation functions in the medium and high z region. The low z behavior of both 
results also fail to catch the feature of the empirical curves. It remains a challenge for further study. 

The indirect ones for kaons are depicted in Fig. m The NLChQM results for zD^ , zD!^ and zD!^ are almost 
identical to zD^ . However, unlike the case of zD^ , NLChQM results of zD!^ , zD^ , and zD^ agree with DSS 
curves excellently between z ~ 0.2 and z — 1. At low z regime NLChQM curves overshoot a little bit but still in 
a reasonable good agreement to the empirical curves. Notice that our results are all within the uncertainty band 
except in the extremely low z region. The last channel we discuss is zDf . Remarkably, the result of zD^ is almost 
identical to the one oi u -^ K. Again our model agrees with DSS curve in the high and medium z region and is above 
the empirical curves in the low z regime. 

V. SUMMARY AND OUTLOOK 

In this article, we have investigated the quark-jet contribution to the fragmentation functions of the pions and the 
kaons using NLChQM and evolute them to Q^=4 GeV^. The current-quark masses, {mu,md,ms) = (5, 5, 150) MeV 
has been used in our calculations. We summarize the important observations in the present work as follows: 



For the direct pion fragmentation functions, NLChQM results agree with the empirical data quite well except 



10 



0.30 
0.25 
0.20 

0.10 

0.05 

0.00 
0.0 




k a 




0.00 
0.0 



0.12 




^ 




0.00 
0.0 



1.0 



(bottom panel, right), 



FIG. 7: The fragmentation functions of zD^ (upper panel, right) and zD^ (upper panel, left), zD^^ 

and zDf (bottom panel, left). The dashed lines denote the results of NJL model. The solid lines represent the results of 

nonlocal chiral quark model. The dotted and dot-dashed lines are the HKNS curves and DSS curves, respectively. HKNS 

and DSS are two empirical parameterizations of the fragmentation functions. The uncertainty band is according to HKNS 

parameterizations. 



at extremely low z regime. 

• In the case of direct kaon fragmentation functions, our result of zD^ is underestimated in the high and 
medium z regime. On the other hand, our result of zDf is overestimated in the high z regime. It is out of the 
uncertainty band of HKNS parametrization, either. Nevertheless, our zDf is still between HKNS and DSS 
parameterizations from medium to low z regime. 

• The most successful channel for NLChQM is m — )■ 7r~ . The agreement between the results of NLChQM and the 
empirical curves is excellent except in the extremely low z regime. 

• For other indirect pion fragmentation functions, such as D1 and D^ , our results are mostly within the 
uncertainty band except in the low z < 0.1 region. Besides our results for those channels are very similar to the 
NJL-jet ones. 

• For the indirect kaon fragmentation functions, NLChQM results agree with one of the empirical curve, DSS 
parametrization, within the regime of 0.2 < z < 1. Furthermore these indirect fragmentation functions lie 
mostly inside the uncertainty bands. 

In summary, we have shown that NLChQM provides an excellent framework to calculate the unpolarized fragmen- 
tation functions. The results agree with the empirical parametrizations quite well in most of the channels. There are 
several directions to improve and extend our current calculations. For example, we have not taken into account the 
axial-current conservation in the present framework, which may become problematic for the nonlocal quark-PS meson 
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interactions [55]. In Ref [15] this contribution has been taken into account to modify the quark distribution functions. 
However, it is not straightforward to include this effect into the calculations of the elementary fragmentation functions. 
This work is in progress. Another issue is to include rj and rj' in the fragmented mesons. Furthermore one should also 
include the vector mesons and baryons in the fragmented hadrons. We expect to continue to study more complicated 
fragmentation functions, such as unpolarized dihadron fragmentation functions, Collins fragmentation functions and 
polarized dihadron fragmentation functions and apply our result to extract the transverse parton distributions of the 
proton. 
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TABLE 1: Flavor factors in Eq. ^. 
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